Radiation effects on the electronic structure of bilayer graphene 
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We report on the effects of laser illumination on the electronic properties of bilayer graphene. By 
using Floquet theory combined with Green's functions we unveil the appeareance of laser-induced 
gaps not only at integer multiples of hfl/2 but also at the Dirac point with features which are 
shown to depend strongly on the laser polarization. Trigonal warping corrections are shown to lead 
to important corrections for radiation in the THz range, reducing the size of the dynamical gaps. 
Furthermore, our analysis of the topological properties at low energies reveals that when irradiated 
with linearly polarized light, ideal bilayer graphene behaves as a trivial insulator, whereas circular 
polarization leads to a non-trivial insulator per valley. 
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I. INTRODUCTION 

Among the many promises s park ed by graphene re- 
search during the last few yearJ^^, graphene optoelec- 
tronics is perhaps one of the brightestPf^. From im- 
proved power conversion of energy harvesting deviceJ^ 
to novel plasmonics propertied^, graphene and re- 
lated materials offer an outstanding playground for te 
study of lig ht-matt er interaction with many potential 
applications^HUnm]. 

Recent studies pointed out the intriguing possibility 
of inducing bandgaps in monolayer graphene by illumi- 
nation with a laser fieldi^Hli] ^j^g peculiar electronic 
structure of graphene and its low dimensionality are cru- 
cial for the occurence of this effect. Further studies have 
predicted observable changes in the conduct anc^^^I^^ 
and optical prop ertied ^ , with a strong dependence on 
laser polarizatioiP^EH, setting off many other interest- 
ing studieJi^lt^Il Moreover, the possibility of controlling 
topological insulators with photocurrenta^, as well as the 
emergence of non t rivial laser - induced topological proper- 
ties and edge state^^^^IllIIIl^ ^jig so-called Floquet topo- 
logical insulatord^MIl^ add more relevance to this area. 

Graphene's thicker cousin, bilayer gra phene (BLG), 
has also shown an enormous potentiaP^, allowing for a 
tunable bandgapP^ as required for the operation of active 
devices. Notwithstanding, the studies mentioned in the 
last paragraph were all centered in monolayer graphene. 
Only in Ref.l^Slj the authors proposed irradiated bilayer 
as a vehicle for inducing a valley polarized current. Here, 
we focus on the electronic and topological properties of 
bilayer graphene illuminated by a laser with frequency 
either in the THz or in the mid-infrared range. In the 
THz range, trigonal warping (TW) corrections are shown 
to induce strong modifications in the theoretical predic- 
tions leading, besides qualitative changes in the spectra, 
to quantitative differences in the laser-induced gaps up 
to a factor of two. 

Moreover, we show that a laser field may also lead to 
polarization-tunable topological properties in BLG rang- 
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FIG. 1: (Color online) (a) Scheme of the considered setup 
where a laser field is applied perpendicular to a graphene bi- 
layer. Panels (b) and (c) show the quasienergy Floquet band 
structure for bulk bilayer without and with trigonal warping 
respectively. These plots are along kx direction {ky — 0), 
solid lines are for circularly polarized light with a frequency 
of 5 THz and an intensity of 0.5mW/ fim^ . The unirradiated 
spectrum is shown with dashed lines. 



ing from a trivial insulator to one with properties akin 
those of a topological insulator. Specifically, we show 
that the low energy properties of BLG illuminated by 
circularly polarized light can be described by a simple 
effective Hamiltonian similar to the one of BLG with a 
bias. Our theoretical analysis shows that although the 
system behaves as a trivial insulator in the presence of 
linearly polarized light, switching the polarization to cir- 
cular transforms it into a non-trivial insulator per valley. 
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II. RESULTS AND DISCUSSION 

A. Floquet theory applied to irradiated bilayer 
graphene. 

The unit cell of bilayer graphene with Bcrnal stacking 
has two inequivalent sites labeled as Al, Bl on the top 
layer and A2, B2 on the bottom layer, they are arranged 
in such a way that atom Bl lies on top of atom A2. Using 
the wave-functions ^' = {ipAi,'4'B2,ipA2,'4'Bi)'^ for the K 
valley and "i/ = (V'B2, V'Ai, V'Bi, "0^2)"^ for the K' valley, 
an effective Hamiltonian for the low energy properties is 
given bji^ 
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where ^ = 1(^1) for valley K{K'), tt — + i py, 
V = {V3/2)ajo/h, V3 ^ (V^/2)a73//i, a = 0.246 nm, 
graphene lattice constant, 70 = 3.16 eV, 71 = 0.39 eV 
and 73 = 0.315 eV. The hopping parameter 73 is respon- 
sible for the trigonal warping effects. 

We apply linear ly/circularly polarized light perpendic- 
ular to the graphene bilayer as shown schematically in 
Fig. [iji. The time-dependent field is introduced through 
the substitution k — > k+eA/h, where the vector potential 
is A{t) = A{cos{rtt), cos{nt + 0)), where (/> = (7r/2}for 
linear (circular) polarization. The Floquet theorenpSlMI 
provides an elegant route to handle this time-periodic 
Hamiltonian {H{t + T) ^ H{t) ^ Ho{k + eA{t)/h), 
where T = 27r/ri), it states that the solutions to the 
time-dependent Schrodinger equation can be written as 



'^K{r,t), where <^„(i) = + T) 



is time-periodic, the Floquet states can be further ex- 

and a 
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panded into a Fourier series [t) =^e^ 
substitution in the Schrodinger equation gives: 
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where = ^ dti/(t)e^("-")"* and e„ is the so- 

called quasi-energy. Simple inspection shows that this 
is an eigenvalue equation analog to the one for time- 
independent systems. There are however two main dif- 
ferences: the role of the Hamiltonian is played by the 
so-called Floquet Hamiltonian Hp = H — ihd/dt; and 
the states belong to an extended Hilbert space which is 
the direct product between the usual Hilbert space and 
the space of time-periodic functions with period T. It is 
straightforward to see that H'^p'™'^ = — nHflSn^m- 

This method has been applied to a variety of systems and 
in particul ar to a c fields such as alternating gate voltages 
in graphen j^^ l '^^ l beyond the adiabatic limit. 

The time-averaged density of states (DOS) gives valu- 
able information on the Floquet spectra in a compact 



form and can be calculated as irfl^EH. To such end we 
compute the Floquet-Green function, defined as Gf = 
(el — Hf)"^, from which the time-averaged DOS is ob- 
tained as DOS{e) = -^/m |rr (Gf(£))o,o}' ^^^ere 
(Gi?)o.o stands for the sub-block of the Floquet-Green 
function corresponding to vanishing Fourier index. 



B. Laser-induced modifications of the Floquet 
spectra 

In the following we will analyze the behavior of the 
quasi-energy spectra and the DOS for various laser in- 
tensities, frequencies and polarizatiori^^J. 

While in monolayer graphene trigonal warping (TW) 
introduces small corrections which become noticeable 
only at high energies (~500 meV), in the case of bilayer 
graphene these corrections are stronger at low energies 
where they lead to a splitting of the Dirac point into a 
structure with four pockets'^^ as shown in the inset of 
Figj2]i. Here we show that these effects, that were ne- 
glected in previous studies of irradiated bilayer, are in- 
deed very important for radiation in the THz range. 

Figure [1] shows the quasi-energy dispersion along a par- 
ticular k direction without (b) and with (c) the TW cor- 
rection in the presence of the electromagnetic field. The 
dashed lines in each figure shows the unirradiated case. 
The field is expected to have the stronger effect at the 
crossing points which, due to the electron-hole symme- 
try, are located at integer multiples of hn/2 above and 
below the Dirac point, as can be seen in Fig. [Tja and 
[1];. The time-dependent perturbation introduces a non- 
vanishing matrix element between the states at those 
crossings, thereby lifting the degeneracies and opening 
the so called dynamical gaps^^ The gap at the charge 
neutrality point is a higher-order effect and will be ana- 
lyzed in more details later. 

Figure [2] shows the DOS for bilayer graphene in the 
presence of either linearly (b) or circularly (c) polarized 
light (5THz) with (solid-line) and without (dashed-line) 
TW. The DOS in the absence of radiation is shown in (a) 
for reference. Although from the discussion before one 
may expect the main corrections to arise only close to the 
Dirac point, Figs. [2]d and[2j; show that they emerge even 
at the dynamical gaps for radiation in the THz range. 

For linearly polarized light, the DOS in the vicinity of 
Ml/ 2 exhibits a depletion area with a linear dispersion 
and a single point of vanishing DOS. This is similar to the 
case of monolayer graphene found in Ref.^^ and is due to 
the fact that the gap depends on relative angle between 
k and the polarization vector, no gap emerges when they 
are parallel. One can also notice that the roughly linear 
dispersion around the dynamical gaps acquire a struc- 
ture with three narrow features on each side when TW 
corrections are included. This is a consequence of the de- 
formation of the iso-energy lines in the k^-ky plane due 
to the TW corrections (see inset of Fig. [2k). 
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FIG. 2: (Color online) a) DOS for bilayer graphene with the 
laser turned ofT. The continous line corresponds to calcula- 
tions performed including trigonal warping (TW) corrections 
and the dashed line without them. Panel a)-inset depicts 
the iso-energy lines for the dispersion of bilayer graphene in 
the absence of radiation, the strong TW distortion is evident. 
Panels b) and c) show the DOS as defined in the text for bi- 
layer graphene in the presence of linearly (b) and circularly 
(c) polarized light (5 THz) with an intensity of 0.5 mSN / iivn? . 
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FIG. 3: DOS as a function of energy for (a) 5THz, (b) lOTHz 
and (c) mid-infrared 30THz radiation, (a) and (b) correspond 
to a laser intensity of 0.5 mW//i m^, (c) is computed for 
10 mW//i m? . Solid (dashed) lines are for circular (linear) 
polarization. Notice the change in the horizontal scale in (c) . 
The structure induced by the TW becomes smoothened as 
the frequency increases. The inset in (b) is a zoom around 
zero energy. 



For circular polarization, two striking observations not 
reported before should be emphasized: a) There is a gap 
opening at zero energy (which also occurs in the absence 
of TW but is much smaller and cannot be distinguished 
in the figure, see inset of Fig. 1-b); and b) the dynam- 
ical gap (which turns out to be linear in the field in- 
tensity as for monolayer graphene) is overestimated by 
almost a factor two when the TW corrections are not 
taken into account. A key ingredient behind these differ- 
ences is again the breaking of the rotational symmetry in 
the kx-ky plane even for low energies. Though the gap at 
the charge neutrality point would require stringent con- 
ditions (being of about 0.3 meV for a laser intensity of 
0.5mW/micronsquare), the physics described here may 
prompt additional research and experiments that may al- 
low directly or indirectly to unveil it. In contrast, the ef- 
fects described at the dynamical gaps {±Ml/2) are much 
stronger and should be observable in low temperature ex- 
periments. Indeed the dynamical gaps are of the order 
of 5K for lOThz radiation at 0.5mW/ and reach larger 
values (up to 30meV, or 350K) for 30THz radiation for 
a power of a few mW/fi^. 

As one moves to higher frequencies, trigonal warping 
effects become less noticeable, though the gaps may be- 
come larger and therefore easier to observe experimen- 



tally. Figure [3] highlights this for three different frequen- 
cies a) 5 THz, b) 10 THz and c) 30 THz (which cor- 
responds to the mid-infrared range) for linearly (dashed 
line) and circularly (solid line) polarized light. 



C. Effective low-energy Hamiltonian description 
and topological considerations 

Though more academic in nature, we now turn to 
an instructive analysis of the low-energy and topologi- 
cal properties of irradiated BLG. Our main fundamental 
question is: are there non-trivial laser-induced topolog- 
ical states to be expected in bilayer graphene? To such 
end, we are interested in obtaining an effective Hamil- 
tonian to describe low energies electronic properties for 
low values of the light intensity in the spirit of Kitagawa 
and coworker^2ll, "VVe will consider only the process when 
one photon is absorbed (emitted) and then re-emitted 
(re- absorbed), in this case applying the continued frac- 
tion method and retaining only the terms of order 0{F^), 
where F = eA/h and in the following h = 1, the effective 
time-independent Hamiltonian can be expressed as: 
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H,ff = Ho + v^M-h ^)V+i + y+iG(+l, n)V^i (3) 

where V±i — _ff for n ~ m ~ ±1 and G{n,ri) — 
t+nn-Ha represents the propagator of a particle with n 
photons. For circularly polarized light, this results in the 
following effective Hamiltonian: 
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where we have assumed 71 S> 3> e. All the terms in 
the diagonal should be multiplied by a factor 77 = ±1 to 
take into account left or right polarization of the light. 
Strikingly, this effective Hamiltonian resembles the 
Hamiltonian of bilayer graphene with a bias, but there 
are some subtle differences that a careful analysis reveals. 
One may argue that laser illumination introduces three 
ingredients: First, it breaks the intra-layer symmetry 
by introducing a te rm similar to Kane-Mele spin-orbit 
term {F^v^/nf^ (if the layers were decoupled the 
system would have a gap solely due to this term); 
second, it breaks the inversion symmetry between the 
two layers (similar to a potential difference between 
layers), an effect which also opens a gap. And third, 
when a graphene-based system with a gap is exposed 
to circularly polarized light an asymmetry between the 
valleys is expected due to the breaking of inversion 
symmetry, an effect similar to opt ical circular dichroism 
for valleys instead of spins .I^ZEII xhe valley degree of 
freedom can be exploited generating valley dependent 
currents as we argue below. 



The gap at fc = is given by 2 x (- 



0, the 



relative importance of these two terms is set by the fre- 
quency f2: for in the THz range the trigonal warping 
term has a leading impact on the gap as previously no- 
ticed in the discussion of Fig. [2j 

To evaluate the topological properties of this effective 
Hamiltonian we reduce the previous 4x4toa2x2 Hamil- 
tonian which describes the effective interaction between 
the non-dimer sites A1-B2. Considering as before 71 ^ e 
the new effective low energy Hamiltonian is given by: 



H = eo 



fc| - s^fc-e -A 



(5) 



where Eq = (73/70)^71 ~ 4 meV, A = vF"^{ff ■ 
C = 1(-1) for valley K{K') , fc± = (fc^ ± iky f/ko, 
ko = 27371 /(•\/3a7Q), a = 0.246 nm and F , fl are now 
dimensionless parameters given in fcg and Sg units respec- 
tively. The parameter s takes values (1,0) to include or 



not the trigonal warping. For values of <C 71 the sec- 
ond term in the expression of A can be neglected and it 
gives a quite simple dependence of the gap with F and 

n' , Gap=2 X Eof^- This expression shows an excellent 
agreement with numerical calculations in the frequency 
range considered to obtain Eq. [4j 

From this effective Hamiltonian it is straightforward to 
calculate the Berry curvature and the Chern numbei^^S!. 
The curvature is given in polar coordinates by: 



n{k,e) 



e?7A(4fc2 - s2 



2(A2 + k^ + k^s^ - 2k^^s cos 361)3/2 ' 



(6) 



and the integration gives an integer non-zero Chern num- 
ber per valley, a quantum valley-Hall state^^. The Chern 
number has opposites values for the two valleys for a 
given handedness of polarization. A valley current will 
be proportional to the Berry curvatur^^. Therefore a 
change in the handedness implies a change in the direc- 
tion of the valley currents, as the sign of the Berry curva- 
ture changes. It provides an effective way to control these 
valley currents. There have be en som e proposals about 
this subject see for instance Ref.'^^'^. On the other hand 
the structure of the Berry curvature reveals the impact 
of the trigonal warping: For low values of A the shape 
of the curvature shows a central dip with a topological 
charge Q=-l and three peaks away from the center and 
separated 120° with Q=l each; in the K' valley we have 
the opposite behavior. This segregation might have an 
impact on the edge currents of a system based on bilayer 
graphene and energies in the Teraherz rangj^. 

A completely different picture is obtained from irradi- 
ating bilayer graphene with linearly polarized light, fol- 
lowing the same procedure as before, one obtains a gap 
at fc = 0, with a peculiar behavior, it does not depend 



explicitly on 17 neither on 73, Gap=2 x 



F^v^ 



The Chern 



number equals zero in every valley, thus the states are 
topologically trivial. 



III. CONCLUSIONS 

In summary, the effects of a laser with frequencies rang- 
ing from THz to the mid-infrared on the electronic struc- 
ture of bilayer graphene are analyzed, highlighting the 
appeareance of laser-induced gaps and their dependence 
with the light polarization as well as the strong influ- 
ence of trigonal warping corrections. For radiation in the 
THz range, trigonal warping in bilayer graphene tends 
to decrease the size of the dynamical gaps at ±.Ml/2^ 
this is very different from the case of monolayer graphene 
where trigonal warping effects are much weakepi^J. Fur- 
thermore, we obtain a time-independent effective Hamil- 
tonian which serves as a starting point for the determina- 
tion of the topological properties of the associated low- 
energy states. We find that while for both polarizations 
there is a small gap at zero energy, their topological origin 
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is different: The Chern number in the presence of Hnearly 
polarized hght equals zero, a trivial insulator, while it is 
a non zero integer, a quantum valley-Hall insulator, when 
the light is circularly polarized. Though more difhcult to 
observe experimentally than the dynamical gaps, further 
work in this direction may open promising prospects for 
exploiting the valley degree of freedom in graphene-based 



structures. 
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